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TECHNICAL NOTE 4154 

APPROXIMATE CALCULATION OF THE COMPRESS IBL E TURBULENT 
BOUNDARY LAYER WITH HEAT TRANSFER AND 
ARBITRARY PRESSURE GRADIENT 
By Eli Reshotko and Maurice Tucker 


SUMMARY 

An approximate method for the calculation of compressible turbulent 
boundary layer with heat transfer and arbitrary pressure gradient is pre- 
sented. The method involves the momentum integral and moment -of -momentum 
equations as siuplified by using Stewartson's transformation. The 
Ludwieg-Tillmann skin-friction relation is used in these equations in a 
form suitable for compressible flow with heat transfer through application 
of the reference enthalpy concept. A tentative extension of Reynolds 
analogy is suggested for estimating heat-transfer effects. 

The method, as applied to insulated surfaces, is quite well founded 
but, for noninsula ted isothermal surfaces, depends on a number of specu- 
lative assumptions . These assumptions are qualitatively proper, and it 
is hoped that they will yield reasonable quantitative results . The 
detailed application of the method for practical calculations is described. 


INTRODUCTION 

In the absence of any rigorous method, various semiempirical pro- 
cedures have been developed for incompressible turbulent -boundary-layer 
calculations. An excellent description of these procedures is given in 
reference 1. In brief, for flows with zero pressure gradient, the 
Karman momentum integral equation is utilized together with an assumed 
boundary-layer -velocity profile, usually a power-law profile, and one of 
several empirical skin-friction relations . For flows with pressure gra- 
dient, an additional or auxiliary equation is required to account for the 
effect of pressure gradient on the boundary-layer -velocity profile. A 
skin-friction relation compatible with the streamwise pressure gradient 
should be employed. 
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Although some procedures use an empirical auxiliary equation, 

Tetervin and Lin (ref. 2 ) have suggested the moment -of -momentum equation 
as an auxiliary equation. The moment-of -momentum equation is obtained 
by multiplying the integrand of the momentum integral equation by a 
distance normal to the surface and then integrating with respect to that 
distance. 

These' calculation methods generally require the simultaneous solution 
of two differential equations, the Karman momentum equation and the 
auxiliary equation. Maskell (ref. 3) has developed a simpler method in 
which the momentum equation is replaced by an empirically determined 
approximation which is directly integrable and thus determines the momentum 
thickness . The profile shape parameter is obtained from an empirical 
auxiliary differential equation. The Ludwieg-Tillmann Bkin-friction 
formula (ref. 4) is used to calculate the skin-friction distribution and 
to determine an approximate separation point for flows with adverse pres- 
sure gradient. 

The present report describes an approximate method for the calculation 
of compressible turbulent boundary layer with heat transfer and pressure 
gradient. The momentum and moment-of -momentum integral equations for com- 
pressible flow are first simplified by using a form of the Stewartson 
transformation (e.g., see ref. 5). The possibility of using the Stewartson 
transformation for turbulent flow was pointed out by Van Le (ref. 6). 

In fact, Englert (ref. 7) and Mager (ref. 8) have respectively applied 
the Stewartson transformation to the Truckehbrodt (ref. 9) and Maskell 
(ref. 3) methods, which enables the calculation of compressible turbulent 
boundary layers over insulated surfaces. In the present report, the 
Ludwieg-Tillmann skin-friction relation is used in a form suitable for 
compressible flow with heat transfer through application of the reference 
enthalpy concept (ref. 10) . An approximation obtained from reference 3 
for the shear-stress distribution through the boundary layer and the 
power -law velocity profile are used to simplify further the moment-of - 
momentum equation. The location of turbulent separation is identified 
as that where the skin friction, when sensibly extrapolated, becomes zero. 
Maskell has shown this procedure to be in good agreement with experiment 
for flows over airfoils where the form factor increases rapidly near 
separation and thus causes a quick fall of skin friction in a small in- 
crement of longitudinal distance. Where separation does not occur so 
rapidly, the proposed method miay yield large errors in separation location. 

A speculative extension of Reynolds analogy is used to estimate 
heat-transfer effects for isothermal surfaces. This extension was sug- 
gested from the results obtained in reference 11 for laminar-boundary- 
layer flows with heat transfer and pressure' gradient. 
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BOUNDARY -LAYER INTEGRAL EQUATIONS 

A brief description, of the essential steps and assumptions involved 
in the proposed calculation method will now be presented. The bulk of 
the development is for a Prandtl number of 1.0. However, in the section 
entitled Heat Transfer, a modification is suggested to include the effect 
of a Prandtl number different from 1.0. All symbols are listed in 
appendix A. 


The transformed momentum and moment-of -momentum integral boundary- 
layer equations, respectively, are: 




U e dX 


HiCHi+l) 2 ^-!) 

2 


1 2 /“Vs 

+ (Hi +“i)"etry 1*0 
L 0 


dY 



( 2 ) 


The development of these equations, which apply to flows with axial sym- 
metry, Is outlined in appendix B. For two-dimensional plane flows, the 

term ^ in equation (l) vanishes. Equations (l) and (2) describe 

the variation of a "transformed" momentum thickness defined as 



( 3 ) 
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and an "incompressible" form factor 



with the transformed longitudinal variable X for a specified external 
velocity distribution U e (x). A distinction is made between the trans- 
formed and incompres sible form factors sincejthe transformed form factor 
for a compressible boundary layer with heat transfer is (obtained from 
appendix C, ref. 12) 



- H i 



(5) 


Physically, Htr is the form factor of a nonlnsulated boundary layer at 
Mach number zero and includes the variation of density with temperature. 
For boundary layers over insulated surfaces, the transformed and incom- 
pressible form factors are identical. 


In order to solve equations (l) and (2) a shin-friction relation 
must first be selected. The Ludwieg-Tillmann relation (ref. 4) has been 
chosen because of its applicability to flows with pressure gradients. 
This relation, however, does not allow the skin-friction coefficient to 
equal zero and, therefore, is in error near separation. This deficiency 
can be overcome in a practical calculation as shown in the section 
entitled Skin Friction, As extended to compressible flow (Bee appendix 
C), the Ludwieg-Tillmann relation in terms of transformed quantities is 


C f = ,-12— . 0.2*6 e- 1 - 5 ™! 
2 Pe“ e 


’MgSoP-tr Y"^ * t e \ 0.732/ t e \0» 268 
, v 0 / Uref/ 


( 6 ) 


This analysis assumes that the fluid Prandtl number is 1.0. The defini- 
tion of reference enthalpy as obtained from reference 10 is 


bref 
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where for Pr = 1.0, hg W = hQ. The temperature t re £ is evaluated at 

the reference enthalpy. Before the solution of equations (l) and (2) is 
possible, some assumptions are still required concerning the shearing 
stress and temperature profiles within the boundary layer. 

For zero pressure gradient, dH^/dX is known to be negative. This 
requires from equation (2) that 



H i 


H-f + 1 


The present analysis uses the following simple relation, valid for 
1.25 < Hi < 1.40: 

£ Z KD ■ l * M stvt ' 8 > 

This relation was obtained by matching with the equivalent numerical 
values of Maskell (ref. 3). Equation (2) appears usable even for flows 
with large pressure gradients inasmuch as the pressure gradient term in 
equation (2) should then far outweigh the friction term (eq. (8)). This 
is somewhat verified by the agreement between experiment and theory in 
reference 12 where the pressure gradient term alone is used to describe 
shock-induced separation. Equation (8) itself is, of course, incorrect 
near separation where ■+ 0. 


The terms in equation (2) involving enthalpy integrals are evaluated 
by drawing an analogy to their values for laminar flow. The results of 
reference 5 show that (from fig. 7 of ref. 5) a zero transformed form 
factor is obtained for a flat plate at absolute zero surface temperature, 
which implies a zero displacement thickness for Mach number zero flow 
with a surface temperature of absolute zero. This is a consequence of 
the similarity between velocity and temperature profiles for the lamina r 
flat plate. Also to be noted from figure 7 of reference 5 is that the 
curves of transformed form factor for insulated and cooled surfaces are 
almost parallel to each other when ■plotted against a pressure gradient 
parameter. For turbulent flow, it is assumed that the qualitative behavior 
is similar to that for laminar flow. This results in the following 
expression for transformed form factor: 


Htr = 




dY 


tr 


i,fp 


■ E t + (^ - l) h : 


0) 
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where Hi fp is the incompressible form factor over a flat plate. Thus, 

f 0 (jl ' l ) “ ■ (§ ' x ) ^ tl0) 

Using equation (10) and manipulating algebraically with the aid of the 

U /Y\^ 

power-law velocity profile — = yield the following identity: 


fofe - l ) * “ - - l ) 


(11) 


By using relations (6), (8), (10) , (ll), and (B4), equations (l) and 
(2) can be rewritten in terms of the physical longitudinal distance and 
Mach number in the following manner: 


d0 tr 9 tr d^e 


dx 

^i = _ 

dx "Mg dx 


0 tr dMe T /h w \ I 

ar sr L 2 + El + ' T 1 ^. 




Hi(H 1 + l) (%-!)' 


q k . x )r 

J| VK) /L H i + 1 2 V H i " l/V H i,fp + 3 /JJ 


- 0.03 Hi(H? - 1) •— 
y tr 


where 


A = 0.123 e 


- 1 , 


56IE i / M e a O 0 tr' 


.-0.268 / 4. \ u. \c 

) (&) fe) 


V 


t e N^0.732/- fce \3.268 


(13) 


(14) 


For Incompressible flow, equations (12) to (14) become 


d^ ©i dU e 0i ar -1.561H-I fU e e ± \~°- 26& . 

"Hr) (15 > 


and 


** - 0.03 Hi (Hi - 1)- 
dx 


/U e 8i\ -0.268 

1 au e 

Hi(Hi + l)(H?-l)“ 

V v 0 ) 

U e dx 

L 2 J 


( 16 ) 
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Equation (15) is the usual momentum integral equation and is identical to 
that of Maskell (ref. 3) except for the present inclusion of axial sym- 
metry. The moment-of-momentum equation (16) differs only slightly from 
that of Mas he 11 which is 


dHi 

dX~ 


f(H,) 


e 1.561Hi cLU e 


u e 0 i 


0.268 


U e 


dX 


g(Hi) 


The first term on the right side of equation (16) has been matched with 
Maskell's values for 1.25 < < 1.40, while the agreement between the 

last term of equation (l6) and Maskell's empirical relation is seen in 
figure 1 to be quite good. Thus, it is felt that the present method is 
in agreement with that of Maskell for incompressible flow, and, therefore 
it is also expected to agree with the method of Mager (ref. 8) for com- 
pressible flow over insulated surfaces. 

The additional assumptions (eqs . (9) to (ll)) made to include non- 
insulated surfaces are of a speculative nature. However, since they are 
qualitatively proper, it is hoped that they will also yield reasonable 
quantitative estimates of the turbulent-boundary-layer characteristics 
over noninsulated surfaces . 


CALCULATION PROCEDURE 

The quantities generally desired in a boundary-layer calculation are 
the boundary-layer integral thicknesses, skin friction, and heat transfer 
In the present procedure. It is first necessary to calculate the 
transformed-momentum- thickness distribution. Anyone interested only in 
estimating heat-transfer rates can proceed directly to the section 
entitled Heat Transfer after first calculating the trans forme d-momentum- 
thickness distribution. 


Transformed Momentum Thickness 

The following form of equation (12) (derived in appendix D) is used 
for the calculation of transformed momentum thickness: 
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The quantities a 0 and v 0 are ^ respectively, the velocity of sound and 
kinematic viscosity evaluated at the stagnation conditions of the local 
external stream. The temperature t re f is evaluated at the reference 

enthalpy h^ef • The local external Mach number Me nnd radius of axi- 
syrametric body R should be known as functions of x in order to per- 
form the calculation indicated in equation (17). For two-dimensional 

1 2155 

flow, the radius R is a constant. Thus, the quantities R 
inside and outside the integral in equation (17) cancel. The flat-plate 
incompressible form factor Hi^fp may, for most purposes, be taken as 
1.3, approximately that of a 1/7 power velocity profile. 

In equation (17) the distance is the starting point of the cal- 

culation. If xp is taken to be the transition location, then the 
transition value of 0-fcr would be obtained from a laminar -boundary-layer 
calculation (by using, for example, the method of ref. 5). For cases w 

where the flow is predominantly turbulent, it is usually sufficient to 
ignore the laminar region and take the origin of boundary-layer development 
at the leading edge (x x = 0). For this latter situation, equation (17) 

becomes 
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to 
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3.268 


Again, for two-dimensional flow, the quantities B 
(20) should he omitted. 


1.2155 


B l. 2155 ^ 

( 20 ) 

in equation 


Transformed Form Factor 


After the transformed momentum thickness is obtained, the quantity 
H-l can be calculated by numerical integration of equation (13): 


<3% 

dx 


1. dM e |H i (H i + l) 2 (H i -l5' 


Me dx 2 

f 1*0 /L H i + 1 



-0.03 Hi(E? - 1) tA- 

tr 


(13) 


where 


-0.268 


0.123 e 


-1.56XH! /' M e a 0 e tr\ ~ ’ ( t e \ 

\ v 0 / Wef / 


0.732 


3.268 


(14) 


For flow over an Insulated surface, equation (13) can be alternately 
written as 


f(Hi) 


f(Hi,l) 



0.060 A 
Hj_ + 1 0tr 


where the function 


f(Hi) 


4 l/(Hi+l) 

(h| - i) 1/2 (z ± + i) 


( 21 ) 


( 22 ) 


is that used in reference 12 and for convenience is plotted in figure 2 
herein. If friction effects are omitted in equation (2l), there results 
the relation applied in reference 12 to shock- induced separation. 
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The initial value of Hi for flow over sharp or pointed bodies may 
be taken as Hf^fp (suggested as 1.3). For blunt bodies, the value of 
Hi at the stagnation point is somewhere between 1.0 and 1.3. A value 
of 1.1 is suggested. Sample calculations in references 3 and 7 show that 
a poor guess of the initial value of Hi is inconsequential, since the 
form factor will tend to reach its proper value in the first few steps 
of calculation. An error greater than 0.1 In the initial value of Hi 
should be avoided. 


Compressible Momentum Thickness 

As shown in appendix C of reference 12, the compressible momentum 

thickness is obtained by using the relation 

* • 

8 = 0 tr (l 4- MfJ 5 (23) 

Skin Friction 

From equation (6), 

C f = 2A^1 + X-ZJ i (24) 

where the quantity A is given in equation (14) . In an adverse pressure 
gradient flow, the point of turbulent separation is obtained by sensible 
extrapolation to zero of the skin-friction coefficients calculated from 
equation (24) . This is shown in the following sketch: 



Since the Ludwieg-Tillmann formula cannot, give Cf = 0, it is probably 
in error for values of Cf close to zero. Maskell (ref. 3) suggests that 
the extrapolation be started where Cf is falling rapidly and be continued 
to Cf = 0 without an inflection point as shown in the sketch. 
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Compressible Form Factor 

As shown in appendix C of reference 12, the compressible form factor 
can be expressed as 


H « =fcr(l + 3L § J: ^ (25) 

where 

=tr - % + - l) O) 

Equation (25) represents well the experimental data of references 13 and 
14 as shown in figure 3. It is also in good agreement with the curve of 
reference 15 for a l/7 power profile. 


Heat Transfer 


In reference 11, it is shown for laminar flow that the ratio of a 
skin-friction parameter to a Reynolds analogy parameter, which appears 
in the relation for heat transfer, is almost invariant with pressure 
gradient. This ratio was then evaluated for zero pressure gradient. For 
laminar flow, this simplification yields good results except for flows 
over heated and moderately cooled surfaces with large favorable pressure 
gradient. It is assumed herein that the corresponding ratio of skin- 
friction parameter to Reynolds analogy parameter for turbulent flow 

e 1 *561Hi/[ Cf /^ st ) (Pr)] behaves similarly) that is, it is essentially 
invariant with pressure gradient and is well approximated by its zero 
pressure gradient value of e“^’^61^i,fp/2(Pr)^/3, The last expression 

makes use of the Colburn flat-plate analogy Cf = (St/2) [(Pr) The 

Reynolds analogy for turbulent flow is then written: 



-2/3 l*56l(Hj_-H^ fp) 
' e 


(26) 


The variation of Stanton number with Prandtl number of (Pr ) -2 / 5 is 
proper for zero pressure gradient. Although for laminar flow the exponent 
of the Prandtl number varies somewhat with pressure gradient, it is felt 
that for turbulent flow the exponent would vary less. In the absence of 
more precise information, the -2/3 exponent will be used with pressure 
gradient. From equations (6) and (26) the expression for Stanton number 
for a compressible turbulent boundary layer is then 
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St = 


0.X23 


0.732 


't£\ ^ /Me^o£tr\ 

■ V V V 0 ) 


0.268 




(27) 


Since the skin friction has been eliminated from equation (27), the heat 
transfer can be directly estimated once the transformed momentum thickness 
is known . 


The heat -transfer rate to the wall is obtained by using the following 
relation: 


■i. 


*>• 

03 

CP 

ro 


q = (St)p e u e (h aw - h*) 


where 


h 


aw 


= h c 


x «• <*)£ x. 


- 1 


wi 


(28) 


(29) 


SPECIAL CASES 

The results of the present method for zero pressure gradient flow 
and 8 tagnat ion-point flow will be described in detail. Both the two- 
dimensional and axisymmetric cases are considered. The value of Hi fp 
for these cases is taken as 1.3. 


Zero Pressure Gradient Flow 


Flat plate . - For flow over an isothermal flat plate (from eq. (20) ), 


1.2155 

e tr 


0.01173 


/t„ \0.732 

tef?) 


X 


/a 0 \° * 2155 

V^o) 




3.268 


(30) 


While the transformed form factor is more strictly described by the 
relation 


dHi 

dx 


-0.03 


HtOq - i)a 


tr 


(31) 



4652 


EfACA TN 4154 


13 


it usually suffices to assume = H^^p = constant. From equations 
(23) and (30), 



and from equations (14), (24), and (30), 



The Stanton number relation from equations (27) and (30) is 


St = 0.043 



, . -2/3 

(Pr)ref 


(34) 


Cone surface . - The thicknesses and skin-friction and heat-transfer 
parameters will be written out in complete form and will also be compared 
to the flat -plate values . With R = x sin 9 C the transformed momentum 
thickness is given from equation (20) by 


9 


1.2155 

tr 


f t \0.732 

0 .01173 ( 7 —— I x 

Vref/ 


/S-.. 0.2155 n once; /+■,-. -.3. 268 

2 . 2155 ^ ^. 2155 ^ 


(35) 


The variation of transformed form factor is more correctly given by equa- 
tion (3l) but may be taken as constant and equal to H^^p. The compres- 
sible momentum thickness is, from equations (23) and (35), 



( 36 ) 
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or , alternately, for the same external Mach number and stagnation condi- 
tions. 


n 

^coue 

e fp 



0.520 


(37) 


The skin-friction relation for a cone is, from equations (14) , (24), and 
(35), 


-0.22 , ^ . 0 .571 .. „ 0.452 


A, , . U.O /1 .v 

-«»£¥) fe © 


(38) 


and 


^£*5222 „ (0.520)- 0 - 26B ^ 1.192 


'f,fp 

From equations (27) and (35), 

>-0.22 , . v 0.571,. .0.452 


(39) 


st . 0.°5l(“L^) (jT^) (£) (Pr )-J (40) 


and 


Ncone = l 

(St)fp Cf^fp 


cone = ^f^cone = 1>l92 


(41) 


These results for the ratios of cone to flat-plate quantities agree 
very well with Gazley's value of 1.176 (ref. 16) and Van Driest's value 
of about 1.15 (ref. 17). 


Stagnation-Point Flow 

The quantities of interest in stagnation-point flow are momentum 
thickness and heat transfer. For stagnation-point flow it is assumed 
that the velocity along the surface varies linearly with the distance 
from the stagnation point and that static 'and total temperatures are 
equal for the region of interest. This treatment of turbulent stagnation 
flows is certainly somewhat fictional. Turbulent flow will be obtained 
only after the transition point and therefore not at the stagnation 
point itself. Also, because of the assumption of equal static and total 
temperatures in the inviscid flow, only a portion of. the subsonic part of 
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the stagnation region can he treated hy using the simple formulas that 
follow. These relations are included, therefore, only as a guide to 
what might be expected from turbulent flow over blunt bodies. 


Two dimensional . - With the previously mentioned assumptions, equa- 
tion (20) yields for momentum thickness: 



while the Stanton number, from equations (27) and (42), is 

0.571 


St = 


0.044[v 0 (B + I)] 0 ’ 22 (£j 




0.22 


-2/3 

(lhr) re f> 


r 0.44 


(43) 


At the stagnation point, the Stanton number is seen to be infinite. It 
is perhaps more instructive to examine the heat-transfer coefficient, 
since the latter is more directly involved in calculating the heat- 
transfer rates: 


*aw ' 


^=(St)p e Ue=0.044p e (Pr)^ 3 [v 0 (B + l)) 0 - 22 ^) 0 ' 571 (^) 0 ' 78 


.0.56 


(44) 


The heat- transfer coefficient and, therefore, the heat -transfer rate at 
a turbulent stagnation point are zero. This is a consequence of using 
the Ludwieg-Tillmann skin-friction relation at zero velocity where it is 
certainly inapplicable. 

Axisymmetr ic . - For momentum thickness with R = x, equation (20) 
becomes 



( 45 ) 
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and the heat-transfer coefficient from equations (27) and (45) is 



-2/3 

0 .044. p e (Pr) re £ [v 


0 (B + 2 .21S5)]°- 22 (^ 0 ’ 571 (^) 0 ‘ 78 


.0.56 


(46) 


These stagnation region results, both two dimensional and axisym- 
metric, are in very close agreement with the results of Van Driest (ref. 
18) . The velocity gradient du e /dx that appears in equations (43) to 
(46) is evaluated. for supersonic flow with M® > 2 from the relation 


du e = fo_ 4 /T 
dx R ap fr 


1 - 


Psp/ 


(47) 


which was derived for bodies having finite^ nose curvature based on the 
assumption of a modified Newtonian flow in the neighborhood of the 
stagnation point. 


CONCLUDING REMARKS 

An approximate method for the calculation of compressible turbulent 
boundary layer with heat transfer and- arbitrary pressure gradient has 
been presented. The method as applied to insulated surfaces is quite 
well founded, but for noninsulated isothermal surfaces it depends on a 
number of speculative assumptions. These assumptions, however, are 
qualitatively proper, and it is hoped that they will also yield reasonable 
quantitative results. The proposed method of estimating the location of 
separation for flow with adverse pressure gradients should be adequate 
for flows where the separation is characterized by a rapid drop of the 
skin friction over a short longitudinal distance. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, September 12, 1957 
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APPENDIX A 


<u 

LO 

to 


A 

a 

B 

Cf 

g(%)) 

H 

h 

M 

«e 

N 

Pr 


SYMBOIS 

skin-friction parameter, eq. (14) 
velocity of sound 

exponent of Mach number in eq. (17) 
local skin-friction coefficient 

functions in Maskell’s moment -of -momentum equation (ref. 3) 

form factor, §*/6 

enthalpy 

Mach number 

local external Mach number 
power of power -law profile 
Prandtl number 


P 

q. 

R 

St 


pressure 

heat -transfer rate 

radius of axisymmetric body 

radius of curvature at stagnation point of blunt body 
Stanton number, 



t temperature 

U transformed longitudinal velocity 

u longitudinal velocity 

V transformed normal velocity 


v 


normal velocity 
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X transformed coordinate along surface 

x coordinate along surface 

Y transformed normal coordinate 

y normal coordinate 

X ratio of specific heats 

A upper limit of integration in momentum of moment -of -momentum 
integral equations 

8* displacement thickness 

9 momentum thickness 

8 C cone half -angle, deg 

v kinematic viscosity 

p mass density 

r shear stress in boundary layer 

Y stream function 
Subscripts : 

aw adiabatic wall 

e local flow outside boundary layer (external) 

fp flat -plate value 

i incompressible 

ref evaluated at reference condition 
s local stagnation value 

sp stagnation-point value 

tr transformed 

w wall or surface value 

0 free -stre am stagnation value 


SS9* 
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1 initial condition 

«> free-stream (static) value 


in 

■J3 

■<* 
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APPENDIX B 


DERIVATION OF TRANSFORMED BOUNDARY-LAYER INTEGRAL EQUATIONS 


For steady compress ib le flow with axial symmetry and a Prandtl number 
of 1.0, the time-averaged continuity and momentum equations may be written. 


Continuity: 


Momentum: 


5(pu) . d(pv) . pu dR 

+ t as 


du , 3u Sp , 

pu y- + pv V- = - y * 1 + y— 

ox ay ox ay 


(Bl) 


(B2) 


The radius of the body R is a function of x alone. For a two- 
dimensional flow, the R term in the continuity equation is not involved. 


Stewartson's Transformation 

The velocities in equations (Bl) and (B2) can he replaced through 
the definition of a stream function: 



(B3) 


The coordinates x and y are transformed by the following form of 
Stewartson's transformation: 


X » 


Y = 



(B4) 


ZSCit 
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The transformed coordinates are now represented by capital letters (X and 
Y ) , and the subscript e refers to local conditions at the outer edge of 
the boundary layer (external). The subscript 0 refers to free-stream 
stagnation values. 


Applying equations (B3) and (B4) to equations (Bl) and (B2) results 
in the following system: 


% + V Y + |I = 0 


UUy + VU v « U e U e _ ^ + 


X Iiq 


Pe /' 

*0 V 


ae 

a o; 




_ St 
2 SY 


(B5) 

(B6) 


where the stream function has been replaced by the transformed velocities 



The relation between physical and transformed longitudinal velocities is 


U = 


u 


fo 

ae 


The integrand of the momentum integral equation is obtained by sub- 
tracting equation (B6) from the product of equation (B5) and the quantity 
U e - U. This results in 


[u(U e - U 


)] x + [v(U e - U)] y + U ex ' 


U e - u + IT e ^ 


•>)] 


(\J - u) — 

R dX 


Pe / a e\ 

P ° PO V®0/ 


__ St 

2 SY 


(B8) 


If equation (B8) is integrated with respect to Y between the limits 
Y = 0 and Y => A where A is a constant distance normal to the trans- 
formed surface sufficiently large so that the conditions U = U e and 
h s = hQ can both be satisfied,- there results 
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d0 t r ®tr d ^e 

ax + u e ax 


2 + Hi + 


rfe - o 

&tr 


dY 


9 tr dR 

r ax " 



where 


*fcr 



( 3 ) 


and 


% 



(4) 


The quantity is the transformed form factor for the boundary layer 

over an insulated surface. For a noninsulated surface, the transformed 
displacement thickness is written: 

' ^■/ A [( 1 '6 + fe' 1 )] aY (B9) 

J 0 

while the transformed momentum thickness is still given by equation (3) 
(see appendix C, ref. 12). 

The moment-of -momentum equation is obtained by multiplying equation 
(BIO) by Y and then integrating with respect to Y between Y * 0 and 
Y « A: 
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(BIO) 


The procedure of reducing this equation to usahle form is that of Tetervin 
and Lin (ref. 2) except for the present consideration of noninsulated 
surfaces. This involves evaluating the integrals in equation (BIO) vxth 

the aid of the power-law assumption U/U e =* (y/a) n and then reverting 
to the form factor % by using the power-law relation % » 1 + 2N. 

The result is 


dX 


1 dU e r + l) 2 ^ - 1) 

Ue dX" L 2 


2 

1 + (% + l)0 t r 



dY 


2(% - 1) 

Hf(H ± + l)8| r 



+ 




(% + 1) 



( 2 ) 


An important result here is that the moment -of -momentum equation is in- 
dependent of radius variation so that equation (2) holds for both two- 
dimensional and axisymmetric flow. This was pointed out in reference 2. 
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APPENDIX C 

ADAPTATION OP LUDWIEG-TILEMANN TURBULENT SKIN-FRICTION 
RELATION TO COMPRESSIBLE FLOW 

The Ludwieg-Tillmann skin-friction relation for incompressible flow 
is from reference 4: 




~ oac -1.561H -a -0.268 
0.246 e Reg 


(Cl) 


where Reg ** UqS/v. The adaptation of this relation to compressible flow 

follows the procedure used by Eckert (ref. 10 ) for the Schultz -Grunow 
relation. With the skin friction expressed in terms of a Reynolds number 
based on x, the fluid properties in the Bkin-friction relation are evalu- 
ated at a reference enthalpy. This reference enthalpy of reference 10 
was selected so as to give good agreement with the skin-friction results 
of many investigations over a wide range of flow conditions. The first 
step then is the conversion of equation (Cl) to x dependence. This will 
be done by using the momentum Integral equation with zero pressure gra- 
dient. The reference enthalpy method is then applied, and the resultant 
relation is reconverted to a Reynolds number based on 6 , again by using 
the momentum Integral equation with zero pressure gradient. The zero 
pressure gradient momentum Integral is 


d 0 C f 

dx 3 2 


0.123 e" 1 - 5613 


vy 0 - 268 


Integrating equation (C2) with constant E yields 


( 02 ) 


0 


r a 

\ -0.268 "1 

(1.268) (0. 123 )e- 1 * 561H 

f) -J 


1/1.268 


Substituting equation (C3) in equation (Cl) then yields 


(03) 



0.364 e -1 * 252 


’u e x\" 0.211 


(04) 
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-3 
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By following the procedure of Eckert (ref. 10), 


T w 


1 2 
2 ^ref 11 © 


-1.23E, 
= 0.364 e ^ 



(C5) 


or 


' c w 


1 

2 


Pe^ 


-1.23B^_ 

0.364 e 



(C6) 


The subscript ref indicates that the quantities should be evaluated at 
the reference enthalpy defined by the relation 

" 5 (sr + x ) + °- 22 (sf - x ) (C7) 

It is customary to demonstrate the adequacy of a compressible skin-friction 
adaptation by comparing the resultant ratio Cf/Cf ± with experimental 
variations. The quantity Cf is the value of skin-friction coefficient 
that would have been obtained without the compressibility adaptation. 

Thus, from equations (C4) and (06), 

Cf ^ / v e y 0 - 211 t e 
Cf,i \ v ref/ ^ref 

For flow over an insulated surface with the adiabatic wall temperature 
equal to stagnation temperature (Prandtl number, 1.0) with viscosity pro- 
portional to temperature and y = 1.4, equation (C8) becomes 

Cf , o,-0.578 , , 

= (1 + 0.144 M|) (C9) 

Cf , i 

This relation agrees very well with the relations of Eckert (ref. 10) and 
Tucker (ref. 15), as seen in figure 4, which in turn are in excellent agree- 
ment with available experimental data. 

For use in the present method, equation (C6), which now includes 
compressibility effects, is reconverted to a Eeynolds number based on 
momentum thickness. The resultant relation is 
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Cf s 


“v 


1 2 
2 Pe«e 


0.246 



(CIO) 


Since , In the present method, the equations are solved in a transformed 
plane (as transformed by the Stewartson transformation) , equation (CIO) 
is transformed similarly with the assumption of a linear viscosity- 
temperature relation. The result is 


Of* 


0.246 e 


2 p e^ 


0.732/, v 0.268 


-1.561Z i /'M e a 0 e tr V°- 268 /t e \ 0 * 732 / ' t 

\w; \t 0 


(6) 


While the assumptions made in the adaptation described here are many, 
it is felt that the result is an adequate description of skin-friction 
coefficient for a turbulent boundary layer with pressure gradient. 


f 
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APPENDIX D 


FURTHER rNTEGRATION' OP MOMENTUM INTEGRAL ELATION 
The momentum integral equation (12) is 


ag tr e tr ^ 
dx + Mg dx 


6 . 


( 2 +Htr ) + 


R dx 


(Dl) 


where 


0.123 e 


-1.56181 J°' 868 ^ *e 


V 


and where 


®tr - El +1^-1) =i,fp 


By following Maskell (ref. 3), let 


® a Ofcr 


and 


M e a O 0 t3 
v 0 


M e a O 0 tr 




J 


(14) 


(9) 


(D2) 


The substitution of equations (D2) into (Dl) yields 

S + ( 1 + k )|S-( 1 + k >[?-i|s i ( H ^ + frl)] <W 

Maskell found that, for incompressible flow, the right side of equation 
(D3) could be represented by 


0 dU e 

0.01173 -4.2^ f 


(D4) 
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with k taken as 0.2155. Since, for incompressible flow, Htr “ H i> i*® 
second term of equation (D4) implies that 

(i + *) (ii + |i|) - 4.2 0 *) 

It then follows that 

(1 + k) ^H ± + (jjj - H i,fp + T+4\ = 4,2 " 1,2155 (h^ " 1 ) H i.fP 

(D6) 

so that equation (D3) may be written as 


g + 1.2155 §g- 0.01173 


/ + \ 0.732/a. \ 3.268 

{&) ft) 


_0_ ^e 
B Mg dx 


(D7) 


where 


B » 4.2 + 1.2155 l ^ 


^w 


*0 


- 1 


(D8) 


The values of B given by equation (D8) with * 3- 3 8X6 in veT V 
good agreement with those for laminar flow (ref. 5) and those estimated 
by Skopets (ref. 19) fo!r turbulent flow. 


Equation (D7) is a first-order linear equation for ® and can be 
integrated to yield 



where 


e a O 0 tr\ 0,2155 


(DIO) 


0 a 0-fcr 
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Figure 2. - Function for calculation of form-factor variation for 
insulated surface. 
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Figure 4. - Comparison of present variation of skin friction with Mach number with 
variations of references 10 and 15 . - 


NACA - Langley Field, Va. 



